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ABSTRACT 


The rotor system of the proposed XH-2/CCR (Circulation 
Control Rotor) prototype aircraft and the state variable 
format of the airframe equations of motion are described. 
Through a study of the eigenvalues and eigenvectors of the 
basic airframe, the effects of uncoupling and cross-coupling 
the helicopter equations of motion were analyzed. The 
control matrix for the CCR was then generated using an 
applicable pneumatic lead angle for Coanda blowing.  Feed- 
back gains, to give the aircraft acceptable flying qualities, 
were calculated for the Stability Augmentation System (SAS) 
and finally the aircraft's sensitivity to changes in pneu- 
matic lead angles were studied. The programs used in the 
above study (capable of handling up to 10x10 matrices) 


are listed for both the IBM-360 and HP-9830 mini-computer. 
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LIST OF SYMBOLS AND ABBREVIATIONS 


Acceleration due to gravity 


Rolling moment about the x-axis due 
to aerodynamic torques 


Pitching moment about the y-axis due 
to aerodynamic torques 


Yawing moment about the z-axis due 
to aerodynamic torques 


Roll rate, angular velocity about the 
x-axis (positive right wing down) 


9 = مص 


Pitch rate, angular velocity about 
y-axis (positive nose up) q = 6 


Yaw rate, angular velocity about 
Z-axis (positive nose right) 


Linear perturbation velocity along 
X-axis (positive forward) 


Linear steady-state velocity along 
the x-axis (positive forward) 


Linear perturbation velocity along 
y-axis (positive out right wing) 


Linear perturbation velocity along 
z-axis (positive down) 


Aerodynamic force along x-axis 
(positive forward) 


Aerodynamic force along y-axis 
(positive out right wing) 


Aerodynamic force along z-axis 
(positive down) 


Phase angle of control system 


Blade aximuth angle from aft position 





1۳ INTRODUCTION 


In the early 1970's the David W. Taylor Naval Ship 
Research and Development Center (DTNSRDC) initiated a 
research program into the feasibility of incorporating a 
Circulation Controlled Rotor (CCR) system ina Navy heli- 
copter. The concept of the CCR and of improving airfoil 
lift-to-drag ratios using Coanda flows had been proven 
earlier by some of the world's leading aerodynamicist's 
and as early as 1959 Dorand had published in the Journal of 
the Helicopter Association of Great Britain an article on 
the application of a jet flap to control a helicopter rotor 
(Ref. 1). Studies and tests continued throughout the 1960's 
with more papers published in both the United States and 
in Europe on the improved performance and possible 
applications of a CCR (Ref. 2). 

The Aviation and Surface Effects Division of DTNSRDC 
continued the research with further tests involving detailed 
pressure measurements of ችንን ኔን ሌዓ | elliptical sections 
in their 15 x20-inch subsonic wind tunnel. These tests 
reconfirmed that extremely high lift-to-drag ratios could 
be achieved by tangentially injecting air through a slot 
in the trailing edge of an airfoil. The results of these 
tests were incorporated in both two-and four- bladed model 
rotor systems for evaluation in the DTNSRDC 8 x10 foot 


windtunnel (Ref. 3). 





The incorporation of a CCR System in a full-size 
aircraft could conceivably offer other advantages over the 
conventional rotor system, in addition to the potentially 
improved aerodynamic performance traits. The conventional 
rotating mechanical swashplate system would be replaced 
by a non-rotating pneumatic plate in a plenum chamber 
located in the blade hub region. Collective control would 
be accomplished by changing the plenum chamber pressure, 
which increases or decreases the Coanda blowing equally at 
all blades via the individual blade supply or collector 
tubes. Cyclic control would be provided by tilting the 
pneumatic swashplate so that there is an azimuthal variation 
in Coanda blowing in each blade. This variation in blowing 
is a result of the changes in volume of air allowed to the 
collector tubes because of changes in the gap between the 
collector tubes and the non-rotating swashplate. This 
non-rotating swashplate and variations in Coanda blowing 
would replace the mechanical cyclic feathering required by 
conventional rotors and therefore eliminate the vibrations 
caused by this one-per-revolution cyclic mechanical move- 
ment. Another important point to recognize is that the 
CCR concept, with the Coanda blowing, will dictate a tor- 
sionally stiff rotor blade or rigid rotor system. This 
is a result of the disparity between the two lift generation 
centroids. The center of pressure due to Coanda blowing 
is near the blade midchord region, while the blade aerodynamic 


center remains near the rotor blade quarter chord point. 


l0 
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The proposed helicopter, with a simple hub and rigid 
rotor blades free of flapping and lag hinges, would result 
in a relatively clean aerodynamic hub system. The reduction 
in rotor and hub drag would be beneficial to the helicopter 
from a performance standpoint. The reduction in moving 
parts in the hub and blade system would also mean a quieter 
helicopter with a lower vibration level than that of a 
conventional rotor system. This latter effect has a favor- 
able potential of improving the "ilities" (maintainability 


and reliability) for helicopter operations. 
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II. BACKGROUND 


Early in 1973 the Naval Air Systems Command (NAVAIR) 
contracted the Lockheed Aircraft Company and Kaman Aero- 
space Corporation to investigate the feasibility of devel- 
oping a full-scale flightworthy Circulation Controlled 
Rotor demonstrator aircraft. In the summer of 1974, some 
twelve months later, both companies returned reports to 
DTNSRDC and NAVAIR stating that: "the concept, while | 
innovative, is completely safe in operation" (Ref. 4) and 
"that there is no fundamental flaw or deficiency in the 
CCR concept and that construction of a full scale CCR 
helicopter is feasible and practical" (Ref. 5). Lockheed 
Aircraft proposed the use of its L286/CCR while Kaman 
suggested "that the Kaman/Navy H-2 aircraft is an ideal test 
vehicle for the CCR concept" (Ref. 4). 

With the additional goal of being able to incorporate 
a CCR system on any off-the-shelf helicopter with no major 
airframe or equipment changes, NAVAIR awarded a contract 
to Kaman Aerospace to "develop, Palia and test" a prototype 
CCR vehicle incorporating the use of the Navy/Kaman H-2 
aircraft. This technology demonstration aircraft will 
tentatively be designated as the Navy XH-2/CCR. 

Preliminary studies by Kaman promoted the belief that 
acceptable flying qualities would be sustained with the 


installed Stability Augmentation System (SAS) of the Kaman 


1 





H-2 with only minor changes in the gains of the feedback 
amplifiers (Ref. 5). Acceptable flying qualities does not 
necessarily mean all stable roots of the aircraft motion 
modes, since a weak oscillatory instability with a time- 
to-double amplitude of greater than three seconds can be 
tolerated by a proficient rotor-wing aircraft pilot. The 
objectives of this research was to confirm that Kaman's 
beliefs were in fact true and to find a suitable feedback 
law for the SAS of the XH-2/CCR such that within the air- 
craft's flight envelope the aircraft flying qualities will 


be acceptable to the evaluation pilot. 
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11 ۳ DISCUSSION 


The study of the helicopter flight dynamics were con- 
ducted using the conventional non-dimensionalized state 
variable format of the aircraft linearized equations of 
motion (Ref. 6), modified to allow coupling of the longi- 
tudinal with the lateral-directional motions. This modi- 
fication is a fairly elementary record-keeping operation 
when using state vector formulations. The basic plant 
matrix, A, of the aircraft linearized equations of motion 
in the state variable format is given in Appendix A. 

The stability derivatives for the XH-2/CCR airframe 
were computer generated by the contractor using the 
MOSTAB-HFA program (Ref. 7) modified for the pertinent 
EN cteristics of the SH-2F airframe and the XH-2/CCR 
main rotor system. The flight conditions analyzed were for 
1.0g level flight at sea level standard conditions. The 
aircraft gross weight was given as 11,000 pounds and a 
rotor tip speed of 615 feet per second (267 RPM) was used 
throughout the calculations. Stability derivatives were 
generated for ው eeds OE: choc hover), 35, 72, 110, 
and 130 knots. These derivatives were computed in May of 
1976 and then updated in November of the same year. The 
calculations made in this research effort are based on 


the updated, November 1976, data. 
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A. THE PLANT MATRIX 

The plant matrix, A, was developed for the longitudinal 
and lateral-directional components and then the fully- 
coupled equations of motion using the contractor generated 
Stability derivatives. The plant matrix, A, was partitioned 


into: 


where [A11] represented the coefficients of the long- 
itudinal stability derivatives and [A551 represented 
the coefficients of the lateral-directional stability 
derivatives. The cross-coupling stability derivatives 
] and d 


were represented by the coefficients of [A, > [A54 


The homogeneous form of the state equations took on 


the conventional form of: 


4 4 è 
H 
۱ 
۱ 
۱ 
| 
ሠ 
4 
ዘ 
ሯ 
| 
۱ 
۱ 


where the state vector, x, included the partitioned 


longitudinal airframe state vector: 


[u, W, Ur 8 ]‏ = ا 
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and the lateral-directional airframe state vector 


respectively: 


x, = IP, r, v, 0] 


Normally in airframe dynamics the fact that most "fixed 
wing" aircraft approach symmetric conditions enables one 
to eliminate the cross coupling terms [A, 5] and [A51] 
and then to represent the aircraft by the decoupled equations 
of motion in the longitudinal and lateral-directional modes. 
Neglecting these cross-coupling terms allows analysis of 
the individual fourth-order systems and the size and shape 
of the applicable longitudinal and lateral-directions modes 
changes only slightly when the more complicated computations 
are made for the fully coupled eighth-order systems. Unfor- 
tunately, helicopters do not enjoy these conditions of 
symmetry and the effects of completely cross-coupling the 
longitudinal and lateral-directional equations are quite 
Significant. 

The fact that the helicopter had large cross-coupling 
effects was shown by the significant change in the eigen- 
values for the uncoupled and fully-coupled computations 
and modal identification could not be completed from the 
results of the cross-coupled forms alone. The conclusion, 
therefore, was that the stability problem could not be 


solved by studying only the uncoupled components of the 


l6 





equations of motion but must be undertaken using the 
fully coupled, eighth-order equations of motion. 

Using the Basic Matrix Control Theory (BASMAT) com- 
puter program (Ref. 8) modified for the HP-9830 mini- 
computer and the IBM-360 digital computer located at the 
Naval Postgraduate School's W. R. Church Computer Center; 
calculations were made at the above stated airspeeds and 
the eigenvalues and eigenvectors of the plant matrix, A, 
were obtained for both the uncoupled and fully-coupled 
equations of motion. The revised programs for both the 
IBM-360 and the HP-9830 are included as part of this paper 
in Appendix B and C respectively. 

The existence of a longitudinal unstable root in the 
uncoupled form was confirmed at all calculated speeds and 
the identification of the different modes in the fully- 
coupled eighth-order system was completed by slowly intro- 
ducing the cross-coupling derivatives into the aircraft 
equations of motion. The eigenvalues and eigenvectors were 
traced and modes identified by letting ور‎ and و2‎ 
equal zero in the eighth-order system and then slowly 
increasing their values until they reached the final values 
of the basic plant matrix, A. The values of the uncoupled 
eigenvalues from the fourth-order solutions and the final 
eighth-order system results are included in Table I. 

A modified root-locus is shown in Figures 1 and 2 
for the airspeed conditions of hover and 130 knots respec- 


tively. The trajectory of the roots is shown as the amount 
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TABLE I 


Eigenvalues of Basic XH-2/CCR Airframe 


Note: Uncoupled airframe listing shows four longitudinal roots first 
followed by four lateral-directional roots. 


Hover (uncoupled) (coupled) 

A= 0.00159 + i 0.10077 A = -0.28206 + i 0.01601 
= -0.23713 = 0.01867 
= -2.64654 = -0.05189 
= 0.02949 + 1 0.27626 = -0.00634 + 1 0.20903 
= -0, 44778 = -5.18781 + i 2.53260 
= -7.71637 

35 Knots (uncoupled) (coupled) 

A = -0.00971 + i 0.16470 A= 0.01046 + i 0.13193 
= -0.30077 = -0.24564 
25537101 = -4,50717 + i 2.34599 
= -0.24771 + i 0.81760 = -0.27651 + i 0.85642 
= -0.05339 = -0.06194 
= -6.61402 

72 Knots (uncoupled) (coupled) 

እ = -0.79763 A = -0.05862 + i 0.15912 
= -0.14599 = -0.27382 
= -0.11719 = -4,40993 + 1 2.33811 
= -2.19846 = -0.42956 + i 1.53113 
= -0.41632 + i 1.42319 = -0.05781 
= -0.05296 
= -6.10178 

110 Knots (uncoupled) (coupled) 

A = -0.38109 + i 0.26126 A = -0.15433 + 1 0.25617 
= 0.31503 = 0.36151 
= -2.92894 22105264955 302.0343] 
= -0.54931 + 1 1.97156 = -0.53360 + i 2.10332 
= -0 09 = -0.04816 
= -5.59397 

130 Knots (uncoupled) (coupled) 

እ = -0.30712 + i 0.29749 À = -0.07151 * i 0.20758 
= 0.31242 = 0.41211 
ع‎ - 22-9 = -4 35708 + 1 2 
= -0.63733 + i 2.22591 = -0.56744 + i 2.34661 
= -0.05636 = -0.08227 
= -5.00383 
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of cross-coupling increased from 0.0 to 1.0 with the latter 
limit corresponding to a fully cross-coupled system. The 
Dutch-roll and long period oscillatory roots are omitted 
for sake of clarity, but their values do not vary 
significantly with cross-coupling as shown in Table I. 

The uncoupled موه‎ tor lateral-directional roots 
may be identified as spiral and roll subsidence roots by 
both the mode shape and time constants using familiar 
analogies from fixed wing aircraft. The remaining two 
non-oscillatory roots from the uncoupled longitudinal 
degrees of freedom would normally correspond to a short 
Period situation in fixed wing aircraft. 

It was observed that the low time constant-real roots, 
one each from the longitudinal and lateral-directional 
degrees of freedom respectively remain almost invariant 
with the amount of cross-coupling until they reach the 
neighborhood of full cross-coupling. At that time the 
longitudinal root becomes weakly unstable with a time-to- 
double amplitude of approximately 37 and 1.68 seconds at 
hover and 130 knots respectively. The latter situation 
definitely required improvement by means of stability 
augmentation. These two real roots could have been expected 
to coalesce into an oscillatory pair as cross-coupling 
Ed. but possibly the close proximity to the almost 
invariant oscillatory long period and Dutch-roll roots 


prevented this action from occurring. 
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The second real-longitudinal root and the roll subsi- 
dence root may be observed to coalesce into a pair of 
complex conjugate (oscillatory roots) almost mid-range in 
the cross-coupling. The coalesced roots have a mode shape 
similar to a pendulum type of motion, but it was noted that 
this mode (in the fully-coupled situations) was quite 
heavily damped. 

The aircraft mode shapes are defined in Table II for 
hover and 130 knots velocities in both the uncoupled and 
fully coupled situations. Without these mode shapes, 
there would be difficulties involved in interpreting the 
Characteristic root migrations as shown in Figures 1l and 
2. The symbol (O) is used on the figures to identify a 
longitudinal plunging subsidence mode which will be spotted 
in Table II with the (w) velocity perturbation being the 
dominant term. This root, which in conventional airframe 
systems would be combined with the rocking mode root (Ll) 
to yield the conventional oscillatory short period mode, 
is the "culprit" which becomes unstable in the fully coupled 
Situation. As will be noted in Figures 1 and 2, the longi- 
tudinal rocking and the lateral roll subsidence modes 


coalesce into an oscillatory pendulum type mode. 


B. THE CONTROL MATRIX 
In addition to the definition of the airframe plant, 
the contractor provided the MOSTAB generated control 


matrices. The derivative information, presumably the result 
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TABLE II 
MODE SHAPE SUMMARY 


A. Hover Conditions 
l. Zero Cross-Coupling 
a. Longitudinal-oscillatory long period 


1» 0.0016 +1 0.1008 Period = 62.3 sec To = 435.94 sec. 
1.0000; arg. 0.0 deg. 
0.0196; arg. -116.8 deg. 


0.0003; arg. 83.8 deg. 
17.1 arg. -5.4 deg. 






OQ SC 


b. Longitudinal-plunging subsidence 


A= -0.2371 T, = 2.92 sec. 
u 0.0700 2 

w {_)} 1.0000 

q 0.0000 

8 0.0000 


C. Longitudinal-fore/aft rocking subsidence 


እ = >2, 647 T, = 0.26 sec. 
تا‎ 1.0000 2 

w \_ J-0.0090 

q -0.1830 

8 0.0690 


d. teral-Directional-Dutch Roll 
X= 0.0295 + i 0.2762 Period = 22.8 sec. T, =) 25.4 sec, 


p 0.0025; arg. 166.4 deg. 
e 0.0250; arg. - 34.4 deg. 
V 1.0000; arg. 0.0 deg. 
d 0.0090; arg. 82.4 deg. 
e. Directional-Spiral subsidence 
እ = -0.448 Ti = 1.55 sec. 
p 0.0044 2 
2109010510) 
V 1.000 
d 0.0099 
f. Lateral-Roll Subsidence 
እ = >7, 716 13 = 0.090 sec. 
p 1.0000 zi 
r( j-0.0294 
V 0.8669 
d 0 





TABLE II (Continued) 


A. 2. Full Cross-coupling 
a. Ilongitudinal-oscillatory long period 


A = -0.2821 + 1 0.0160 Period = 392.7 sec. 


1.0000; arg. 0.0 deg.) T1 7 2:46 sec 
0.7102; arg. 14.1 deg. 2 

0.0025; arg. 173.4 deg. 

O 0059, arg. = 3.3 deg. 
0.0008; arg. 176.0 deg. 
0.0490; arg. 174.5 deg. 
0.1924; arg. 182.6 deg. 
1 aro. =. 0.7 deg. 










ዘ 


DOSE‏ ۳۵ < همه 


Si 


Cross-coupled Pendulum Mode 
ሙንን) ያቄ ን 6 Period = 2.48 sec. 


og. 1 Sec‏ 188 مھت 
arg. 155.40 deg.‏ ;0.0798 
arg. 2.14 deg.‏ ;0.4219 
arg. -155.84 deg.‏ ;0.0731 
arg. - 42.92 deg.‏ ;0.7618 
arg. 177.95 deg.‏ ;0.0386 
arg. 0.0 deg.‏ ;1.0000 
arg. 163.10 deg.‏ ;0.1320 







NJ 


e 4 ۲۷۲۲ عم و مرو‎ ው 
1 


Lateral-Direction - Coupled Dutch Roll Mode 
x = -0.0063 + 1 0.2090 Period = 30.1 sec. 


5 










u) f1.0000; arg. 0.00 deg.) 71 = 107-26 Sec. 
۷ 0.0162; arg. -171.87 deg. 2 
q 0.0014; arg. 3.04 deg. 
9 (_/0.0065; arg. - 88.70 deg. 
p 0.0006; arg. 145.88 deg. 
6 0.0126; arg. 1.69 deg. 
V 0.4099; arg. 33.20 deg. 
0.0028; arg. 122.38 deg. 
d. Lateral- 
= -0.0519 T. = 13.36 sec. 
1.0000 2 
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I3. 


TABLE II (Continued) 


e. Longitudinal 
à = 0.0187 


e HÜ pu £ نم‎ 





130 Knots 


Ir 


Zero Cross-coupling 


2 


Tj, = 37.13 sec. 


a. Longitudinal-oscillatory Long Period Mode 
De 0880718 12022975 Period = 21.1 sec. T. = 2.26 sec. 


u 075812; arg. -67/.70 deg. 
W 1.0000; arg. 0.0 deg. 
q 0.0026; arg. 39.45 deg. 
0 0.0061; arg. -96.46 deg. 






b. Longitudinal - 


A= 0.3124 

u 1.0000 
w(_ 20.7721 
q -0.0043 
8 0.0158 


C. Longitudinal - 


እ = -3.0256 

ኒጊ 0.0185 
wt{_) 1.0000 
q -0.0108 
ሀ 0.0036 


d. Lateral-Directional Dutch Roll Mode 


> 


0.0185; arg. 0.02 deg. 
1.0000; arg. 0.0 deg. 
0.0107; arg. 179.96 deg. 
0.0036; arg. = 0.0 deg. 






e d Hr 
| 


29 


12 


T 


Dol be 


= 063735 1 2.2259 Peraed = 2.8 sec. T 


1 


NI 


= 2.22 sec. 


= 0.23 sec. 


= 1.09 sec. 


N| 





B. 


2: 


TABLE II (Continued) 


e. Lateral- 
à = -0.0564 
p -0.0077 
er. 0.0203 
V 1.0000 
b 0.1363 
f. Lateral- 
A = -5.0038 
p 1.0000 
ri ۱-1 
v{ ۱-64 
ፅ -0.1999 


Ti = 12.29 sec. 
2 

> 0.14 sec. 
2 


Full Cross-coupling 


a. Longitudinal-oscillatory Long Period Mode 
እ = -0.0715 + 10.2076 Period = 30.3 sec. T, = 9.69 sec. 


© > HO od SC 


1.0000; 
0.4987; 
0.0014; 
0.0064; 
0.0045; 
0.0027; 
0.3260; 
0.2057: 


مو إنم 


arg. 0.0 deg. 
arg. 10.05 deg. 


arg. 34.18 deg. 
arg. - 74.983 deg. 
arg. 167.65 deg. 
arg. 81.22 deg. 
arg. 25.33 deg. 
arg. 883.34 deg. 


b. Lateral-Directional - Coupled Dutch Roll Mode 
À = -0.5674 + 1 2.3466 Period = 2.7 sec. T, =1.22 sec. 


CO SN GO mo = ር. 


ር 


(00.27 
0.2045; 
0.0023; 
0.0009; 
0.0045; 
0.0110; 
1.0000; 
0.0019; 


l 
arg. -109.36 deg. p 
arg. ፦ 20.66 deg. 
arg. - 74.39 deg. 
arg. 98.99 deg. 
arg: 159.81 GEG. 
arg. - 76.806 deg. 
arg. 0.0 deg. 
arg. 56.22 deg. 
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TABLE II (Continued) 


c.  Cross-couple Pendulum Mode 
À = -4.3571 x i 2.8196 Period = 2.2 sec. T, = 0.16 sec. 






il 
0 0.0177; arg. 12.86 deg. 2 
W 1.0000; arg. 0.0 deg. 
q 0.0214; arg.  -66.01 deg. 
ፀ8ሏሊ .7 0.0041; arg. - 3.95 deg. 
pí )0.0412; arg. 82.02 deg. 
Y 0.0066; arg. - 8.08 deg. 
V 0.2695; arg. 37.04 deg. 
ፅ 0.0079; arg. -65.07 deg. 
d. Longitudinal- d. Lateral 
A= 0,4121 T, = 1.68 sec. A = =0;0823 T4 = 8,42 sec. 
u -0.7732 u 1.0000) 2 
wW 1.0000 W 0.0781 
q 0.0006 q -0.0000 
2 0.0014 6 أ‎ _([ 9 
p -0.0031 pf )|-0.0041 
r -0.0014 ገ 0,0076 
ሇ -0.0342 V 0.5428 
d -0.0075 4 04 
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of pneumo-dynamic modeling in the hub and blade blowing 
sections, was supplied as airframe forces and moments per 
unit (3060 psf) pressure variation at the pneumatic swash- 
plate referenced to the aircraft axis. The cyclical varia- 
tion of the plenum pressure upon control was given by: 
C (2) .. The coefficient of the cosine ytype 
variation in the control matrix. 
C (3) .. The coefficient of the sine y type 
variation in the control matrix. 
where, U, represents the blade azimuthal angle. One can 
visualize the longitudinal and lateral cyclic controls as 
rotating the pneumatic swashplate about an orthogonal set 
of axes that leads the blade azimuthal angle by some angle, 
$, in an analogous manner to the practise on mechanical 
swashplates. Then, as shown in Figure 3, the apparent 
longitudinal and lateral cyclic control matrices become, 


by a coordinate rotation, as follows: 


۳۰ 1) 


Geco C (3) Sin ò 


)2( ظط 


۳ ۲ Et ‘cos 4 


where B(1) and B(2) are the nal and lateral cyclic 
control matrices respectively. No attempt was made to relate 
these control matrices to actual control stick motions, 
although reasonable estimates could have been made. 

The criteria employed in selecting the control lead 


angle was the following: 
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FIGURE 3.  VISUALIZATION OF THE PLENUM 
PRESSURE LEAD ANGLE 


AIRFRAME LONGITUDINAL 
BODY AXIS 
+X 


E EXE Eve 
E e AXIS 


Va DET دا‎ ) ۳2 2۳۸۱ LATERAL 
BODY AXIS 
ንም ምጭ en 
He ት 


لذ 


Gs, età CEG 
cc . ር CONTROL AXIS 





9ھ 


. کس‎ AZIMUTHAL ANGLE) 


SWASHPLATE 


LEAD ANGLE) — CYCLIC CONTROL LAW 
TRANSFORMATION - 


5፡1 | کی‎ 4 pd fce 
B(2) m cos SEN 
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O Longitudinal cyclic .. produced negligible 
rolling moment. 


O Lateral cyclic .. produced negligible 
pitching moment. 

Although absolute satisfaction of these constraints 
concurrently with one choice of lead angle, è, was not 
physically realizable, it was remarkable that a single 
value of lead angle, $ = 40 degrees, provided a satisfactory 
solution in an engineering sense. 

Table III lists the variation of pneumatic lead angle 
versus airspeed for satisfaction of the longitudinal and 
lateral cyclic constraints respectively. The selection of 
forty (40) degrees as an engineering answer is in accord 
with independent analysis done by Kaman Aerospace. The 
tabulation of the 8 )1( and B(2) control matrices for a 
pneumatic lead angle of forty (40) degrees are presented in 
Table IV. inspection of the fifth row of B(1) and third 
row of B(2) provides an indication of the reasonableness 
of the solution. 

The effects of plus and minus five (5) degree changes 
in pneumatic lead angle will be described during the analy- 
sis. This type of sensitivity analysis will provide an 
indication of the airframe stability root sensitivity in 


the compensated mode. 


C. THE FEEDBACK LAW 
The study of the impact of various feedback control 


laws for reducing or removing the longitudinal instability 
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TABLE III 


Calculations of the Pneumatic Lead Angle vs. Airspeed 





B(1) = C(2) cos 0 = C(3) sin ۵ 
C(2), 
Bl = arctan SI ey 
B(2) = C(2) sin ¢ = C(3) cos d 
-C(3)q 
95% = arctan EOF 

HOVER 951 = 42,322 

de = 439250? 

= 6 

35 KNOTS و‎ 3716 
$52 - mE 

72 KNOTS 951 = 40,872 

E 6 

"53 33.09 

110 KNOTS 93 = 39.83° 

- 0 

Ge 36563 

130 KNOTS 951 = 40.612 

= : o 

Ge 42.56 
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TABLE IV 


Coefficients of the Control Matrix Bix 


and B for pneumodynamic lead angle 


2x 
of (40) degrees 


HOVER 35 KTS 72 KTS 1 ከ ከ 130 KTS 















— 1220916 211 213 -8.474 FU 5993 SILL 













- 13 9 295195 22 2 =27561 





o 








6.726 





35356 78565 51 n 229 














0.0 0.0 0.0 0.0 





0.0 








0.174 





08295 20:205 





1.164 1124 














=05055 -1.446 0.905 በ S 






20:267 















2.842 3.049 2.030 





2.816 





4.897 





0.0 


SI 5.486 3.860 2.086 0.639 1.711 
22 -0,765 SaL. 072956 11.21 11.53 
E32 20572 0.084 0.837 0.425 0.304 
B42 0.0 0.0 0.0 0.0 0.0 
E52 28.72 ھ7‎ ۹۹ ۹۹7 18.30 MED 
B2 25 CU c 0:521 0.965 0.585 
E12 12.94 10.36 9.244 10.60 9.996 
E32 0.0 0.0 0.0 0.0 0.0 
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began with an attempt to vary the available longitudinal 
feedback gains and to investigate the eigenvalues and 
eigenvectors as these gains were varied. 

The modified plant matrix, A', was developed in the 
traditional matrix manner in the uncoupled and fully- 


coupled state variable format where: 


A' = A - Bk 


When only longitudinal cyclic control is considered, the 
control effectiveness matrix [B] becomes an eight-by-one 
matrix while the feedback gain coefficient matrix [k] 

becomes a one-by-eight matrix. The matrix product, Bk , 


is an eight-by-eight matrix. 


B Ou, EER, V 0] 


k = Ik. 7 K 7 Ko ; SÉ ۰ KL ; Ku ; Kal 
The earlier confirmation that the unstable root was primarily 
| associated with the longitudinal airframe modes was the 
reason for only employing feedback in the longitudinal 
cyclic control when developing the modified plant matrix, 
[A']. An arbitrary set of moderately damped oscillatory 


stable roots were selected: 


À => IOS 2 1 3.21 
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which corresponded to: 


undamped natural frequency = 8.38 sec - 


E 
li 


dimensionless damping ratio - 0.2 


ሆኘ 
Il 


Since it had been established in the basic plant matrix 

that the short period mode was the dominant instability, 

the arbitrarily selected second-order system (a form of 
modal control) was applied to the augmented two-by-two 
matrix, A', and the closed form solution calculated to yield 
values of k and KG: Gain values of Ky = 0.05 and zm = 0.22 
were determined to yield the desired results but the applica- 
tion of these gains alone, to the uncoupled four-by-four 
matrix and the fully-coupled eight-by-eight state variable 
problem failed to produce favorable results and the insta- 
bility remained with the aircraft matrix SE all of the 
calculated speeds. 

The search for the acceptable feedback law using gains 
of Ky and SE continued using the HR-9830 with a further 
modified BASMAT program that would automatically search for 
acceptable feedback gains in the range of K and Kg equal 
to minus one (-1) to plus one (+1.0). Although values of 
E and Ka could be found that would drive the augmented 
matrix stable (negative real parts of the eigenvalues) 
at each calculated speed, these values were not sequentially 
related to speed and, furthermore, were random in nature, 
often changing signs more than once as speed increased from 


hover. 





The decision was then made to employ pitch attitude 
and pitch rate (ko and Kg) feedback respectively based on 
the knowledge that both pitch and pitch rate information 
were presently available in the H-2 aircraft. Programs 
existed for both the IBM-360 and the HP-9830 for accom- 
plishing this search for acceptable feedback gains of k 


q 


and k but while the IBM-360 was much faster in the actual 


gr 
computations (approximately twenty seconds of CPU time were 
required for the computations of the eigenvalues and eigen- 
vectors for one set of plant, control and feedback gain 
coefficients vice twenty-five minutes for the HP-9830), 
the HP-9830 allowed for a much more convenient search. The 
HP-9830 allowed the programmer to make "in-line" decisions 
on changes in the feedback gains based on the previous 
results with approximately thirty minutes between output 
results. This removed the problem of the long delays 
| encountered because of the turnaround time of the IBM-360. 
Typical turnaround times were two to five hours depending 
upon the computer usage at the time of program input. This 
long turnaround time was a result of the low job priority 
assigned the program by the computer center, a. result of the 
complexity of the program and the large amount of core 
memory required. 

The rationale of feeding back Ka and k, proved fruitful 


6 


for values of: 


Ka = 0.45 and ka = 0.85 
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or the feedback matrix taking the form of: 
Fa رز را‎ 0.85, 0, 0, 0, 01 


producing stable eigenvalues at all speeds calculated with 
the exception of hover. These gains did leave an unstable 
oscillatory root in hover, with a time to double amplitude 
of 18.15 seconds, over six times the three second minimum 


time to double amplitude allowed. 


D. PNEUMATIC LEAD ANGLE SENSITIVITY CHECK 

Concern had been shown over the possibility that 
variations in the pneumatic lead angle (5) blowing could 
cause dramatic changes in the stability characteristics of 
the aircraft. It was recognized that changes in pneumatic 
lead angle would result in direct changes to the control 
matrix, B. In order to study the sensitivity to changes 
in lead angle, new values of the control matrix were com- 
puted for pneumatic lead angles of thirty-five (35) and 
forty-five (45) degrees. The results of these calculations 
are listed in Table VI. The IBM-360 program was further 
modified to compute the new eigenvalues and eigenvectors 
using the basic plant matrix, A, the computed feedback gains, 
k, given N and allowing the longitudinal pneumatic 
control matrix, B, to vary for values of ọọ equal to 


thirty-five and forty-five degrees. The computed results 
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TABLE V 


Eigenvalues of the augmented matrix, A', at 
velocities of: 


Hover, 35, 72, 110 and 130 Knots 


Hover: 

= -0.00019 

= -0.24473 

= -0.47775 

= -0.97459 n 

= 0.03819 - i 0.27676 

= -6.78748 + i 3.60274 
35 Knots: 

= -0.01709 

= -0.04164 

= -0.24762 + i 0.81412 

= -0.67857 + i 0.13139 

= -5.67357 + i 3.42617 
72 Krots: 

= -0.01445 

= -0.07094 

= -0.63074 + i 0.37790 

= -0.43457 +. i 1.47940 

= -5.45309 + i 3.07803 
110 Krots: 

= -0.00326 

= -0.04764 

= -0.49037 + i 0.55734 

= -0.54836 + 1 2.06899 

= -5,61868 + i 3.46986 
130 Knots: 

= -0.01401 

= -0.07035 

= -0.39270 + i 0.61107 

= -0.59046 + i 2.34511 

>ታ፡.ፒ...]፣፤፣ን፣ TT ሪህ ጋ 
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showed that the variation of the pneumatic lead angle by 
plus or minus five (5) degrees had very little effect on 
the stability characteristics of the augmented plant 


matrix. 
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IV. CONCLUSIONS 


A study has been made of the basic stability traits of 
the Kaman Aerospace Corporation XH-2/CCR helicopter, which 
is presently being constructed under NAVAIR contract as a 
technology demonstrator for the Circulation Control Rotor 
concept. 

The airframe was defined by contractor generated 
stability and control derivatives which were then used to 
develop eigenvalues and eigenvectors for the system. The 
plant matrix (which characterizes the airframe) was generated 
by the MOSTAB program as modified to accommodate the CCR 
system, and the matrix coefficients represent the output 
from the program when it was operating in the 18 degree-of- 
freedom situation, i.e., six airframe degrees-of-freedom 
plus flapping, torsion and lead/lag degrees-of-freedom for 
the four blades (6+3x4 = 18). The uncompensated airframe 
characteristic roots at hover are in close accord with 
results obtained by the contractor for the six degree-of- 
freedom airframe, including a mild non-oscillatory insta- 
bility (t, 7 37.1 sec.) that has been identified as due to a 
longitudinal short period type mode. 

A feedback law has been identified using pitch attitude 
and pitch rate feedback into the longitudinal cyclic control 
that provides reasonable characteristic roots for the 


airframe.  Presumably, further improvements could be obtained 
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by providing feedback in the lateral cyclic control system. 
The effect of varying the cyclic control lead angle on the 
pneumatic swashplate was investigated and found to be slight. 

It should be noted that a unique feedback control law 
is not possible in modal control theory, when multiple 
Control inputs (longitudinal and lateral cyclic) are 
Available. Another way of stating this fact is that it is 
possible with several control inputs (and feedback laws) to 
have the same eigenvalues, but with different eigenvectors. 

Finally, the characterization of the eigenvectors and 
identification of the eigenvalues with relevant modes was 
made possible by using a form of root locus analysis where 
the prime parameter was the amount of cross coupling. The 
trajectory of the characteristic rcots as cross-coupling 
was linearly altered provided a physical insight into the 
history of the various roots. 

Future studies are suggested to include estimating the 
mechanical gearing to the cockpit controls and then obtaining 
airframe response time histories for selected control inputs 
such as stick doublet type motions. Time histories can be 
generated quite readily using principles from control theory 
combined with calculated relations for system's transition 
Matrices. It is quite possible that time history calcula- 
tions for the compensated airframe will provide a better 
guide for selecting the control law. In any event, such 
Studies will enhance satisfaction as to the question of 


airframe response behavior being reasonable. 
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APPENDIX B 


IBM-360 MODIFIED "BASMAT" PROGRAM 
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